USAARL  Report  No.  95-13 


Design  of  Digital  Low-pass  Filters 
for  Time-Domain  Recursive  Filtering 
of  Impact  Acceleration  Signals 


By 

Nabih  Alem 

Aircrew  Protection  Division 
and 


Matthew  Perry 
University  of  South  Alabama 


19950412  072 

January  1995 


Approved  for  public  release;  distribution  unlimited. 


United  States  Army  Aeromedical  Research  Laboratory 
Fort  Rucker,  Alabama  36362-0577 


Notice 


Qualified  requesters 

Qualified  requesters  may  obtain  copies  from  the  Defense  Technical  Information  Center  (DTIC), 
Cameron  Station,  Alexandria,  Virginia  22314.  Orders  will  be  expedited  if  placed  through  the 
librarian  or  other  person  designated  to  request  documents  from  DTIC. 

Change  of  address 

Organizations  receiving  reports  from  the  U.S.  Army  Aeromedical  Research  Laboratory  on 
automatic  mailing  lists  should  confirm  correct  address  when  corresponding  about  laboratory 
reports. 

Disposition 

Destroy  this  document  when  it  is  no  longer  needed.  Do  not  return  it  to  the  originator. 

Disclaimer 

The  views,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should 
not  be  construed  as  an  official  Department  of  the  Army  position,  policy,  or  decision,  unless  so 
designated  by  other  official  documentation.  Citation  of  trade  names  in  this  report  does  not 
constitute  an  official  Department  of  the  Army  endorsement  or  approval  of  the  use  of  such 
commercial  items. 

Reviewed: 


ibviN  T.  MASON 
LTC,  MC,  MFS 
Director,  Aircrew  Protection 
Division 


Released  for  publication: 


DENNIS  F.  SHANAHAN 


Colonel,  MC,  MFS 
Commanding 


Unclassified 

:E<iURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


la.  REPORT  SECURITY  CLASSIFICATION 
Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


lb.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 
Approved  for  public  release,  distribution 
unlimited 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S)  5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

USAARL  Report  No.  95-13 


6a.  NAME  OP  PERFORMING  ORGANIZATION  1 6b.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 
U.S.  Army  Aeromedical  ResearcH  Of  •PPhcblt)  ^  Medical  Research  and  Mat 


Laborator 


6c  ADDRESS  {City,  State,  end  ZIP  Code) 
P.O.  Box  620577 
Fort  Rucker,  AL  36362-0577 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


MCMR-UAD-CI 


U.S.  Army  Medical  Research  and  Materiel 
Command _ 


7b.  ADDRESS  (C/fy,  State,  end  ZIP  Code) 
Fort  Detrick 

Frederick,  MD  21702-5012 


8b.  OFFICE  SYMBOL  I  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  eppliceble)  I 


8c  ADDRESS  (City,  Stete,  end  ZIP  Code) 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PI 

ELEMENT  NO.  I  N 


e2787A 


IWORK  UNIT 
■ACCESSION  NO. 


87A818  ED 


1 1 .  TITLE  (IrKlude  Security  Cleuificetion) 

Design  of  digital  low-pass  filters  for  time-domain  recursive  filtering  of  impact 
acceleration  signals 


12.  PERSONAL  AUTHOR(S) 

Nabih  M.  Alem  and  Matthew  Perry 


13a.  TYPE  OF  REPORT  13b.  TIME  COVERED  14.  DATE  OF  REPORT  (Year,  Month,  Oayl  15.  PAGE  COUNT 

Final _  FROM _ TO _  1995  January  30 


16.  SUPPLEMENTARY  NOTATION 


COSATI  COOES 


GROUP 


01 


02 


18.  SUBJECT  TERMS  {Continue  on  reverse  if  necessary  end  identify  by  block  number) 
SUB-GROUP  I  digital  signal  processing,  recursive  digital  filters, 
impact  testing 


19.  ABSTRACT  (Continue  on  reverse  if  necessery  end  identify  by  block  number) 

The  purpose  of  this  report  is  to  present  simplified  and  portable  computer  pro¬ 
grams  to  design  and  implement  recursive  digital  Butterworth  filters.  The  design 
method  described  in  the  report  uses  the  bilinear  transfoirmation  to  convert  the 
well  established  analog  design  to  digital  formulation.  Fortran  codes  for  the 
design  and  implementation  of  the  filter  are  included  in  the  report.  The  result¬ 
ing  design  is  compared  to  another  filter  proposed  in  SAE  J211  guideline  for 
filtering  impact  acceleration  signals.  Although  both  design  methods  produced 
filters  which  are  within  corridors  prescribed  in  the  J211  guideline,  the  standard 
design  produces  a  Butterworth- like  response  whereas  the  J211  design  does  not. 
Additionally,  the  standard  design  methods  and  the  implementation  algorithms  allow 
the  design  of  any  even-order  filter  which  may  be  required  for  filtering  other 
long-duration  signals. 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
13  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT. 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 
Chief,  Science  Support  Center 


OD  Form  1473,  JUN  86  i 


21.  ABSTRACT  SECURITY  CLASSIFICATION 
□  DTIC  USERS  Unc lassified 


22b  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL 
205-255-6907  MCMR-UAX-SS 


Previous  editions  ere  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

Unclassified 


Acknowledgments 


This  work  was  partially  supported  by  funding  from  U.S.  Army  Night  Vision  and 
Electronic  Sensors  Directorate,  Fort  Belvoir,  Virginia.  The  diligent  review  of  the  manuscript  by 
Mary  Gramling  and  the  meticulous  editing  by  Udo  Volker  Nowak  are,  as  always,  deeply 
appreciated. 


Table  of  contents 


Page 


List  of  figures  .  2 

Introduction  . 3 

Objectives  .  6 

Methods  .  6 

Results  .  9 

Discussion  .  10 

Conclusions  .  11 

Notes  on  fi'equency  response  plots  . 12 

References  .  13 

Appendix  A.  Program  to  filter  in  the  frequency  domain .  22 

Appendix  B.  Program  to  design  Butterworth  low-pass  filters  .  23 

Appendix  C.  Program  to  implement  a  second-order  filter  in  the  time  domain .  26 

Appendix  D.  Program  to  design  a  filter  per  SAE  J21 1  (draft)  guidelines .  28 


,  Isccfssion  lop _ 

PsiT  S  GRA&I  Ef 

DTIC  T.&B  □ 

Unxuinoynced  □ 

J'iis  t  i.1’  i  e  at.  Ion - - — .. 


List  of  figures 


Figure  Page 

1 .  Comparison  of  causal  and  noncausal  time  responses  to  a  step  input  signal  .  5 

2.  Frequency  response  of  100-Hz  filter  designed  for  variable  sampling  rates 

by  standard  Butterworth  method  (top)  compared  to  SAE  J21 1 

CFC  60  filter  (bottom) .  14 

3 .  Frequency  responses  of  300-Hz  filter  for  variable  sampling  rates 

designed  by  standard  Butterworth  method  (top)  compared 

to  SAE  J211  CFC  180  filter  (bottom) .  15 

4.  Frequency  responses  of  1000-Hz  filter  for  variable  sampling  rates 

designed  by  standard  Butterworth  method  (top)  compared 

to  SAE  J2 1 1  CFC  600  filter  (bottom)  .  16 

5.  Frequency  responses  of  1650-Hz  filter  for  variable  sampling  rates 

designed  by  standard  Butterworth  method  (top)  compared 

to  SAE  J2 1 1  CFC  1 000  filter  (bottom) .  17 

6.  Frequency  responses  of  1 00-Hz  filter  for  1 0  kHz  fixed  sampling  rate 

designed  by  standard  Butterworth  method  (top)  compared 

to  S  AE  J2 1 1  CFC  60  filter  (bottom) .  18 

7.  Frequency  responses  of  300-Hz  filter  for  10  kHz  fixed  sampling  rate 

designed  by  standard  Butterworth  method  (top)  compared 

to  SAE  J21 1  CFC  180  filter  (bottom) .  19 

8.  Frequency  responses  of  1000-Hz  filter  for  10  kHz  fixed  sampling  rate 

designed  by  standard  Butterworth  method  (top)  compared 

to  SAE  J21 1  CFC  600  filter  (bottom) . 20 

9.  Frequency  responses  of  1650-Hz  filter  for  10  kHz  fixed  sampling  rate 

designed  by  standard  Butterworth  method  (top)  compared 

to  SAE  J21 1  CFC  1000  filter  (bottom) . 21 


2 


Introduction 


The  U.S.  Army  Aeromedical  Research  Laboratory  (USAARL)  often  is  tasked  with  the 
assessment  of  injury  potential  from  impacts  and  jolts.  This  requires  the  analysis  of  accelerations 
and  forces  obtained  from  transducers  mounted  in  human-like  manikins  and  test  forms,  and  generated 
during  impact  tests.  With  the  exception  of  the  signal  conditioning,  the  analysis  is  conducted  at 
USAARL  almost  entirely  on  personal  computers  (PC)  to  perform  the  analog-to-digital  conversion 
of  conditioned  signals,  filter  the  digital  signals,  extract  injury  parameters  from  the  signals,  and 
display  the  results  in  graphical  forms  on  the  PC  screen  and  on  an  attached  printer. 

One  of  the  standards  for  processing  impact  signals  is  the  Society  of  Automotive  Engineers 
J21 1  guidelines  for  instrumentation  for  impact  testing  (SAE,  1994).  The  J21 1  requires  signals  from 
impact  tests  to  be  filtered  using  one  of  four  channel  frequency  classes  (CFC)  of  low-pass  filters  and 
specifies  acceptable  frequency  response  for  each  filter  class.  The  four  filters  are  designated  as  CFC 
60,  1 80,  600,  and  1000.  It  is  clear  from  the  J21 1  filter  specifications  that  they  were  derived  from 
analog  Butterworth  filters  whose  comer  frequency  is  equal  to  the  CFC  designation  divided  by  0.6. 
The  comer  of  a  low-pass  Butterworth  filter  is  defined  as  the  frequency  at  which  the  signal  loses  one- 
half  of  its  power,  i.e.,  where  the  signal  magnitude  attenuation  is  equal  to  v^,  or  -3  decibels  (dB). 


Thus,  the  comer  of  CFC  60  filter  is  at  100  Hz,  CFC  180  at  300  Hz,  CFC  600  at  1000  Hz,  and 
that  of  CFC  1000  at  1650  Hz.  In  previous  versions  of  the  J21 1,  acceptable  roll-off  slopes  of  filters 
ranged  from  12  to  24  dB/octave,  i.e.,  filters  with  2,  3,  or  4  poles  were  acceptable.  The  1994  draft 
proposes  upper  and  lower  slopes  which  are  24  dB/octave,  suggesting  that  a  4-pole  filter  is  the  basis 
for  the  requirement.  The  method  of  filtering  is  left  up  to  the  user  and  may  be  done  with  analog 
filters  or,  as  is  the  current  practice  in  many  testing  facilities,  with  digital  filters. 

Filtering  is,  perhaps,  the  most  critical  phase  in  the  processing  of  impact  signals.  Its  primary 
function  is  to  eliminate  undesired  high-frequency  noise  that  obscures  the  underlying  signature  in  the 
signal.  The  importance  of  filtering  becomes  evident  when  considering  that  filtering  reduces  the 
peaks  in  the  signal  and  peaks  often  are  used  for  assessment  of  protective  devices.  The  proliferation 
of  personal  computers  has  promoted  the  conversion  of  analog  signals  to  digital  ones,  and  increased 
the  need  for  sophisticated  digital  signal  processing  algorithms  to  replace  the  functions  traditionally 
reserved  for  analog  electronic  systems. 

Because  filter  design  formulas  are  well-established  in  the  continuous-time  world  of  electrical 
engineering,  they  often  are  adapted  for  digital  filtering.  Analog  Butterworth  filters  have  the  property 
of  having  a  maximally  flat  frequency  response  in  the  pass-band,  and  an  asymptotic  roll-off  beyond 
the  comer  frequency.  The  roll-off  slope  is  a  function  of  the  order  of  the  filter;  however,  regardless 
of  the  order,  the  attenuation  always  is  -3  dB  at  the  comer  frequency.  More  detailed  description  of 
the  characteristics  of  these  and  other  filters  may  be  found  in  many  textbooks,  e.g.,  Oppenheim  and 
Schafer,  1975. 
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A  digital  filtering  method  which  has  been  used  at  USAARL  is  to  transform  the  digital  signal 
to  the  frequency  domain,  using  fast  Fourier  transforms  (FFT),  then  attenuate  each  frequency 
component  by  an  amount  equal  to  the  Butterworth  function  at  that  frequency.  Since  both  the  real 
and  imaginary  portions  of  the  frequency  magnitude  are  attenuated  by  the  same  amount,  no  phase 
distortion  is  introduced  and  the  resulting  filter  is  phaseless.  Appendix  A  is  a  listing  of  a  Fortran 
subroutine  that  implements  this  filtering  method.  Although  this  method  has  proven  effective  for 
most  applications,  it  has  two  main  disadvantages.  First,  because  the  filtering  is  performed  in  the 
fi’equency  domain,  there  are  restrictions  placed  by  the  FFT  algorithm  on  the  number  of  samples  in 
the  signal.  Thus,  the  signal  may  not  exceed  a  predetermined  length  and,  often,  the  number  of 
samples  must  be  a  power  of  2.  This  means  that  if  a  longer  duration  signal  is  needed  for  some 
analysis,  the  sampling  rate  must  be  reduced  in  order  to  meet  the  limited  size  and  longer  duration 
requirements.  Reduction  of  sampling  rate  may  be  tolerated  up  to  a  point  below  which  events 
containing  high  frequencies  would  not  be  captured  in  the  digitized  signals.  An  example  of  this 
situation  is  the  repeated  jolts  signal  where  several  sharp  impacts  occur  separated  by  time  lapses  that 
increase  the  overall  duration  of  the  entire  signal. 

Second,  and  a  more  serious  disadvantage  of  frequency  domain  filtering  has  to  do  with  the 
causality  of  the  filter.  In  analog  filters,  the  output  signal  is  produced  only  as  a  result  of  an  input 
signal.  It  is  clear  that  output  does  not  "anticipate"  the  oncoming  step,  but  slowly  rises  as  a  result  of 
it.  The  delayed  response  of  this  causal  filter,  shown  in  Figure  1,  distorts  the  phase  relationships 
between  different  signals  and  must  be  removed  in  order  to  synchronize  the  timing  of  events  recorded 
in  various  chaimels.  Phaseless  filtering  is  achieved  in  frequency-domain  FFT  filtering  which 
eliminates  the  time  delay  between  input  and  output.  Another  approach  to  phaseless  filtering  in  the 
time  domain  is  to  filter  the  signal  once  in  the  forward  direction,  then  a  second  time  in  the  reverse 
direction.  It  may  be  seen  fi-om  Figure  1  that  such  noncausal  filter  unfortunately  produces  an  output 
that  anticipates  the  event  and  starts  responding  to  it  before  it  occurs.  The  disadvantage  of  this 
behavior  is  the  distortion  of  preimpact  state  which  is  essential  in  some  applications  where  the  pre¬ 
impact  value  is  used  as  the  "zero  state"  of  the  transducer  output. 

Time-domain  recursive  filtering  addresses  the  disadvantages  FFT  filtering.  First,  with 
recursive  filtering,  we  do  not  have  to  contend  with  FFT  algorithms  that  restrict  the  number  of 
samples  and  sampling  rate.  The  only  limit  is  the  amount  of  memory  which  may  be  set  aside  in  the 
PC  hardware.  More  important,  filtering  may  be  done  only  in  the  forward  direction  to  produce  causal 
filters.  The  user  continues  to  have  the  option  to  produce  a  phaseless  filter  at  the  cost  of  losing  its 
causality.  With  this  flexibility,  time-domain  filters  offer  an  attractive  alternative  to  frequency- 
domain  ones.  This  was  recognized  in  the  newly  proposed  J21 1  instrumentation  guidelines  which 
now  include  an  appendix  that  provides  the  implementation  of  a  phaseless  fourth  order  Butterworth 
filter  (SAE,  1994). 
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input  step  signal 

causal  response  of  analog  filter 

response  of  noncausal  digital  filter 


Figure  1.  Comparison  of  causal  and  noncausal  time  responses  to  a  step  input  signal. 
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Objectives 


(1)  To  develop  a  computer  program  to  design  a  digital  Butterworth  filter  of  arbitrary 
order  and  comer  frequency  for  time-domain  implementation. 

(2)  To  develop  a  computer  program  to  implement  the  digital  Butterworth  filter 
recursively  in  the  time  domain. 

(3)  To  compare  the  frequency  response  of  the  provided  design  and  implementation  with 
that  of  the  proposed  J21 1  filters. 


Methods 


The  method  for  designing  the  desired  digital  filter  is  to  transform  the  known  analog  filter 
equations  into  the  digital  domain  using  the  bilinear  transformation.  Details  of  the  procedure  are 
described  in  many  digital  signal  processing  textbooks  (e.g.,  Cappellini,  Constantinides,  and  Emiliani, 
1978)  and  are  summarized  here.  The  squared  magnitude  function  of  analog  Butterworth  filters  is 
defined  in  the  complex  s-plane  by 


H(s)H(-s)  =  - ^ - 

1  +  (sY 


(1) 


where  N  is  the  order  of  the  filter  and  s  is  a  complex  variable.  These  filters  have  their  poles  in  the  s- 
plane  equally  spaced  on  a  circle  of  radius  equal  to  the  comer  (-3  dB)  frequency, 

The  bilinear  transformation  which  defined  by 


s 


k 


1  -  z 
1  +  z* 


(2) 


is  a  simple  algebraic  substitution  which  is  applied  to  the  Butterworth  filter  of  equation  (1).  This 
yields  the  transfer  fimction 


1 


/  N 

2 

,  1  -  z-1 

+ 

— 

k 

1  1  - 

(3) 
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which  can  be  expressed  in  the  frequency  domain  by  letting  z'^  —  e-^'^  to  give  the  squared  magnitude 
of  the  frequency  response  function 


1 

1  + 

Ptan^- 

2 

N 

(4) 


For  a  Butterworth  filter,  the  squared  magnitude  is  equal  to  one  half  at  the  comer  frequency  regardless 
of  the  value  of  N.  In  particular,  for  N  =  1 ,  the  denominator  of  equation  (4)  must  be  equal  to  2  at  the 
comer  co  =  o,. ,  i.e.. 


I  +  tan^  — -  -  2  (5) 

2 

This  defines  the  constant  k  of  the  bilinear  transformation  for  a  Butterworth  filter: 

=  1  I 


which  may  substituted  in  equation  (4)  to  produce 


\H{en\^ 


1 


1  + 


tan(fa)/2) 

tan(o)^/2) 


(7) 


Equation  (7)  is  the  squared  magnitude  of  the  digital  Butterworth  filter  which  was  obtained 
by  applying  the  bilinear  transformation  to  the  analog  function.  The  real  and  imaginary  parts  of  the 
poles  of  Butterworth  filter,  whose  order  iV  is  even,  are  written  as  follows  (Gold  and  Rader,  1969): 


1  2  2 
m  •'in 


U  = 
m 


(1  - 


2  + 

✓  m 


and 


(1  - 


(8) 


1 


where 


2m +1  ,  w  2m+l 

=  tan —  cos - n  and  v  =  tan —  sm - ■k 

2  IN  2  IN 


(9) 


for  m  =  0, 1, 2, 2A^-1.  Time-domain  implementation  is  done  simply  by  cascading  these  sections 
to  achieve  the  desired  overall  filter. 

The  next  problem  to  address  is  the  method  of  implementation  of  the  filter  in  the  time  domain. 
In  a  cascade  implementation,  filtering  will  be  done  in  stages,  where  output  of  one  stage  is  used  as 
input  to  the  next  one.  Let  X(n)  be  the  n-th  input  sample  in  the  sequence  of  unfiltered  digital  signal, 
and  Y(n)  the  corresponding  sample  in  the  filtered  output  signal.  Then,  the  output  sample  is  given 
by  the  difference  equation: 


Y{n)  =  a^Xin)  +  a^X{n-\)  +  a^X(ri-l)  -  b^Y{nA)  -  b^Y{n-2)  (10) 


where  Qq,  Qj,  bj,  b2  are  the  coefficients  of  a  second-order  filter  sections,  derived  from  the  real 
and  imaginary  parts  of  the  poles. 


The  coefficients  of  the  J21 1  filters  are  essentially  those  of  a  Butterworth  design,  except  that 
the  comer  frequency  is  defined  in  terms  of  the  J21 1  channel  filter  class,  and  an  empirical  factor  is 
introduced  into  the  equations.  As  with  our  Butterworth  design,  the  4th  order  filter  is  achieved  by 
cascading  two  second-order  sections  which,  in  the  J21 1  guideline,  are  identical.  Given  a  signal 
sampled  at  intervals  of  T  seconds  (inverse  of  sampling  rate  in  Hz),  the  five  coefficients  of  a  J21 1 
filter,  whose  channel  filter  class  designation  is  C,  are  given  by: 


= 


(1  +  v^co^  +  wj) 


a^  =  2 


^2  ^0 


6.  = 


6.  = 


2(1-0)^) 


(1  +  v/2a)^  +  (jil) 

(1  -  v/2a)^  +  a>l) 

(1  +  v/2o)^  +  to^) 


(11) 


8 


where 


(12) 


The  (C/0.6)  is  the  comer  frequency  of  the  Butterworth  filter  and  accommodates  the  common 
usage  of  "channel  filter  class"  designation  C  instead  of  a  comer  fi'equency.  The  other  1 .25  constant 
is  an  empirical  constant  which  will  be  discussed  later. 

Finally,  we  will  describe  briefly  the  method  used  for  generating  the  frequency  response 
curves  for  this  report.  The  magnitude  response  at  a  given  fi-equency  may  be  generated  by  passing 
a  sine  wave  of  that  fi-equency  through  the  filter  and  simply  recording  the  amount  of  attenuation 
caused  by  the  filter.  By  judicious  selection  of  fi-equencies,  a  curve  may  be  generated  by  connecting 
all  the  fi-equency  response  points  of  individual  sine  waves.  This  procedure  was  the  basis  for  a 
computer  program  which  was  written  in  Microsoft  Fortran  to  design  filters,  and  to  generate  and  plot 
frequency  response  curves  on  a  personal  computer.  The  program  invokes  the  two  filter  design  and 
implementation  routines  listed  in  Appendixes  B,  C,  and  D.  Although  the  program  allows  the 
selection  of  the  filtering  direction  (forward  or  backward),  only  forward  filtering  was  selected.  This 
was  done  after  demonstrating  that  the  direction  of  filtering  only  affects  the  time  response  of  the 
signal,  i.e.,  time  delay  of  events  in  the  signal,  but  not  its  frequency  response.  The  curves  were 
plotted  graphically  on  the  PC  monitor  display.  Then  a  screen  image  capture  utility  was  used  to  save 
the  plotted  response  to  a  bitmap  file.  Later,  a  Microsoft  Windows-based  utility  (Paint)  was  used  to 
retrieve  each  screen  image  file  and  print  it  on  a  laser  printer. 


Results 

The  design  formulas  given  by  equation  (8)  were  coded  into  the  Fortran  subroutine  listed  in 
Appendix  B.  The  recursive  implementation  algorithm  given  by  equation  (10)  also  was  coded  in  the 
Fortran  subroutine  listed  in  Appendix  C.  The  J21 1  filter  design  formulas  given  in  equations  (11)  and 
(12)  were  coded  in  the  subroutine  listed  in  Appendix  D.  These  three  subroutines  were  used  in  a 
program  ( not  included  in  this  report)  which  was  written  to  generate  frequency  response  plots  for 
user-specified  parameters. 

To  test  the  accuracy  of  the  design  routines,  four  filters  were  designed.  These  are  the  four 
SAE  J21 1  filters  commonly  used  in  processing  anthropomorphic  manikin  transducer  signals  (e.g., 
accelerations,  forces,  moments  ...)  and  obtained  during  impact  and  crash  testing.  The  same  filters 
were  designed  twice:  first  using  the  bilinear  transformation  (Appendix  B),  then  using  the  proposed 
J21 1  formulas  (Appendix  D).  The  resulting  frequency  response  plots  are  shown  in  Figures  2  through 
5.  These  plots  were  generated  under  "ideal"  sampling  rates,  i.e.,  such  that  the  sampling  rate  of  each 
sine  wave  signal  was  at  least  10  times  the  frequency  of  the  sine,  or  at  least  10  times  the  comer 
frequency  of  the  filter,  whichever  was  higher. 
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Ordinarily„however,  an  analog  test  signal  is  sampled  at  a  fixed  rate  which  is  supposed  to  be 
at  least  twice  the  highest  frequency  contained  in  the  signal,  but  usually  is  5-10  times  that  frequency. 
In  order  to  simulate  this  realistic  condition,  a  fixed  sampling  rate  of  10,000  samples  per  second  (10 
kHz)  was  used  to  generate  frequency  responses  of  the  four  J21 1  filters  using  the  two  methods  of 
Appendixes  B  and  D.  The  10  kHz  is  the  sampling  rate  recommended  in  the  J21 1  for  sampling 
impact  test  signals.  Results  of  this  simulation  produced  the  eight  frequency  responses  are  shown  in 
Figures  6,  7,  8,  and  9. 


Discussion 


It  is  clear  from  the  top  graphs  in  Figures  2,  3,  4,  and  5  that  the  standard  design  routine 
(Appendix  B)  and  the  recursive  time-domain  implementation  algorithm  (Appendix  C)  produce  the 
well-known  Butterworth  response  when  the  sampling  rate  is  allowed  to  vary  to  accommodate  the 
frequencies  of  the  sine  waves  being  filtered.  A  standard  Butterworth  filter  has  a  -3  dB  attenuation 
at  the  comer  frequency,  and  an  asymptotic  roll-off  at  the  rate  of  12  dB/octave  for  each  second-order 
section.  In  the  standard  design,  the  asymptote  crosses  the  frequency  axis  exactly  at  the  comer 
frequency,  as  shown  in  the  top  graphs  of  Figures  2,  3,  4,  and  5.  On  the  other  hand,  the  bottom 
portions  of  the  same  figures  demonstrate  that  the  frequency  responses  produeed  by  the  J21 1  design 
formulas  (Appendix  D)  do  not  result  in  the  standard  Butterworth  response,  even  though  they  remain 
within  the  specified  J21 1  response  corridor. 

This  deviation  is  attributed  to  the  1.25  empirical  constant  of  equation  (12).  Recall  that  J21 1 
proposes  to  cascade  two  identical  second-order  sections  to  produce  the  desired  fourth  order  filter. 
However,  by  using  identical  sections,  the  attenuation  at  the  design  comer  frequency  is  no  longer  -3 
dB,  as  expected  in  a  Butterworth  design,  but  doubles  to  -6  dB  at  the  comer  frequency.  Since  the 
frequency  response  is  a  continuously  decreasing  fiinction,  there  exists  a  frequency  where  the  fourth- 
order  attenuation  response  crosses  the  -3  dB  level.  It  is  this  frequency  that  the  constant  1 .25  tries 
to  capture.  By  designing  a  second-order  filter  with  a  comer  frequency  1 .25  times  higher  than  the 
desired  comer  of  the  fourth-order  filter,  the  overall  effect  will  be  to  produce  an  attenuation  of  -3  dB 
at  the  comer  frequency  of  the  desired  filter.  This  is  evident  in  Figures  2  through  5  (bottom  graphs) 
where  the  J2 1 1  design  produced  the  desired  attenuation  at  the  comer.  Unfortunately,  the  intersection 
of  the  roll-off  asymptote  with  the  frequency  axis  does  not  move  to  the  -3  dB  frequency,  but  remains 
at  the  original  comer  of  the  second-order  filter,  where  the  attenuation  is  now  -6  dB.  In  other  words, 
the  flat  portion  of  the  pass  band  does  not  extend  as  far  as  the  standard  design  at  the  new  comer,  but 
starts  rolling  much  earlier. 

Although  an  ideal  Butterworth  response  may  be  achieved  when  the  sampling  rate  is 
unrestricted,  in  reality,  the  sampling  rate  is  limited  by  hardware  and  software  considerations  to  a 
fixed  rate.  For  example,  the  J21 1  guideline  recommends  a  sampling  rate  of  10  kHz.  Using  this 
sampling  rate  to  illustrate  its  effects  on  the  frequency  response,  it  is  clear  the  frequency  response  of 
the  filters  deviate  noticeably  from  the  ideal  Butterworth  response  regardless  of  the  filter  design 
method  (Figures  6,  7,  8,  and  9).  Since  this  deviation  is  unavoidable  when  the  sampling  rate  is  fixed. 
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it  is  necessary  tq  set  limits  for  the  deviations  beyond  which  the  response  would  be  unacceptable. 
The  J21 1  corridors  which  are  superimposed  on  all  the  figures  in  this  report  provide  these  limits.  Of 
course,  these  limits  are  intended  for  processing  force  and  acceleration  signals  from  manikin  crash 
tests,  and  may  be  redefined  when  other  applications  emerge.  In  addition,  the  deviation  occurs  at 
frequencies  two  to  three  times  the  comer  frequency.  In  general,  these  frequencies  should  have  been 
reduced  already  by  the  use  of  antialias  analog  filters  prior  to  the  analog-to-digital  conversion. 

Finally,  the  recursive  filtering  routine  provided  in  Appendix  C,  which  implements  the 
difference  equation  (10),  should  be  discussed  briefly.  Because  two  prior  samples  (n-1  and  n-2)  are 
required  to  compute  each  current  (n)  sample,  it  is  clear  the  first  filtered  sample  that  can  be  computed 
is  point  no.  3.  Therefore,  a  starting  method  has  to  be  devised  to  deal  with  the  initial  conditions  of 
the  filter.  In  the  code  provided  in  Appendix  C,  the  method  used  is  to  extend  the  starting  segment  of 
the  signal  by  reflecting  points  2  and  3  symmetrically  about  point  1 .  This  provides  two  additional 
starting  points  which  are  used  to  start  the  algorithm,  then  discarded.  Other  methods  may  have  to  be 
devised  by  the  user  to  deal  with  or  take  advantage  of  specific  initial  conditions.  Alternatively,  the 
user  may  start  digitizing  the  analog  signal  earlier  than  the  event  of  interest,  then  discard  the  startup 
segment  of  the  filtered  signal  in  order  to  avoid  the  initial  effects  of  the  filtering  process.  Another 
comment  on  the  method  given  in  Appendix  C  has  to  do  with  memory  allocation.  Since  one  of  the 
objectives  of  filtering  in  the  time  domain  is  to  increase  utilization  of  computer  memory  for  long 
signals  instead  of  auxiliary  storage  required  by  FFT  filtering,  this  was  accomplished  in  the  provided 
code  at  a  small  cost  in  the  algorithm  complexity. 


Conclusions 


A  computer  program  was  developed  to  design  a  Butterworth  low-pass  digital  filter  using  the 
bilinear  transformation.  A  companion  program  was  written  to  implement  the  filter  recursively  in 
the  time  domain.  The  code  for  these  two  programs  is  highly  portable  and  may  be  recoded  in  any 
high-level  language.  The  two  programs  offer  a  flexible  and  memory  efficient  alternative  to  FFT 
filtering  and  have  been  demonstrated  to  be  stable  and  accurate.  The  frequency  response  of  the  J21 1 
filters  was  compared  to  those  designed  by  our  methods.  No  advantage  was  found  in  one  method 
over  the  other  when  the  sampling  rate  was  fixed  to  10  kHz.  However,  the  ideal  Butterworth  filter 
can  be  achieved  precisely  with  the  bilinear  transformation  design  program  offered  in  this  report. 
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Notes  on  frequency  response  plots 


•  The  frequency  response  plots  included  in  this  report  were  generated  point  by  point  by 
passing  sine  wave  signals  of  different  frequencies  through  the  filter. 

•  For  the  first  eight  responses  (Figures  2, 3, 4,  and  5),  the  sampling  rate  was  variable,  that  is 
each  sine  wave  signal  was  sampled  at  a  rate  at  least  1 0  times  the  comer  frequency  of  the  filter 
or  at  least  10  times  the  frequency  of  the  sine  wave  being  filtered,  whichever  was  greater. 
For  the  remaining  eight  frequency  responses  (Figures  6,  7,  8,  and  9),  the  sampling  rate  was 
fixed  at  10  kHz  which  is  the  rate  recommended  by  the  J21 1  guideline.  However,  for  the 
purpose  of  generating  these  plots,  the  sampling  rate  was  never  allowed  to  be  less  than  10 
times  the  comer  frequency  in  order  to  allow  a  sufficent  number  of  samples  per  period. 

•  To  deal  with  the  end  effects  of  the  filter,  only  the  middle  third  portions  of  the  input  and 
output  time  signals  were  scanned  to  determine  the  peak-to-peak  span  of  the  sine  wave 
signals. 

•  The  erratic  behavior  in  some  of  the  frequency  response  plots  at  high  frequencies  of  both 
designs  may  be  explained  by  the  small  numerical  value  of  the  peak-to-peak  range  in  the 
filtered  output  which  tends  to  be  overcome  by  numerical  rounding  and  truncation  errors 
resulting  in  the  observed  behavior.  In  an  actual  digital  signal,  these  high  frequencies  should 
have  been  attenuated  already  by  antialias  filters  so  that  these  numerical  artifacts  should  not 
be  significant. 

•  Both  second-order  filtering  stages  were  done  in  the  forward  direction.  Filtering  direction 
affects  the  time-response,  but  not  the  frequency  response  of  the  filter. 

•  The  straight  lines  above  and  below  the  frequency  response  curve  are  those  defined  in  the 
J21 1  as  boundaries  of  the  accepted  frequency  response  corridor.  The  faint  straight  line  in  the 
middle  of  the  corridor  is  the  asymptote  to  standard  Butterworth  response  and  has  a  slope  of 
24  dB/octave,  i.e,  12  dB/octave  for  each  of  the  two  second-order  filter  sections. 

•  The  unlabeled  horizontal  grid  line  between  0  and  -1 0  dB  corresponds  to  the  -3  dB  attenuation 
level. 
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Figure  2.  Frequency  response  of  100-Hz  filter  designed  for  variable  sampling  rates  by  standard 
Butterworth  method  (top)  compared  to  SAE  J21 1  CFC  60  filter  (bottom). 


Figure  3.  Frequency  response  of  300-Hz  filter  for  variable  sampling  rates  designed  by  standard 
Butterworth  method  (top)  compared  to  SAE  J21 1  CFG  180  filter  (bottom). 


Butteruorth  corner  1000  Hz,  stages:  FF 


Variable  sampling 


Figure  4.  Frequency  response  of  1000-Hz  filter  for  variable  sampling  rates  designed  by  standard 
Butterworth  method  (top)  compared  to  SAE  J21 1  CFG  600  filter  (bottom). 


Butteruorth  corner  1650  Hz,  stages:  FF 


Uar table  sanpling 


Figure  5.  Frequency  response  of  1650-Hz  filter  for  variable  sampling  rates  designed  by  standard 
Butterworth  method  (top)  compared  to  SAE  J21 1  CFC  1000  filter  (bottom). 


Butteruorth  corner  100  Hz,  stages:  FF 


10000  Hz  sampling 


Figure  6.  Frequency  response  of  100-Hz  filter  for  10  kHz  fixed  sampling  rate  designed  by  standard 
Butterworth  method  (top)  compared  to  SAE  J21 1  CFG  60  filter  (bottom) 


Butteruorth  corner  300  Hz,  stages:  FF  10000  Hz  sanpling 


0  500  100 

Frequency  (Hz) 


J211  CFC  180  corner  300  Hz,  stages;  FF  10000  Hz  sanpling 


Frequency  (Hz) 


Figure  7.  Frequency  response  of  300-Hz  filter  for  10  kHz  fixed  sampling  rate  designed  by  standard 
Butterworth  method  (top)  compared  to  SAE  J21 1  CFC  180  filter  (bottom). 


Figure  8.  Frequency  response  of  1000-Hz  filter  for  10  kHz  fixed  sampling  rate  designed  by  standard 
Butterworth  method  (top)  compared  to  SAE  J21 1  CFC  600  filter  (bottom). 


Frequency  (Hz) 


Figure  9.  Frequency  response  of  1650-Hz  filter  for  10  kHz  fixed  sampling  rate  designed  by  standard 
Butterworth  method  (top)  compared  to  SAE  J21 1  CFC  1000  filter  (bottom). 


Program  to  filter  in  the  frequency  domain. 


II II  It  II II II  II  II  It  II II II II  II  II  II  II II II II II II  It  It  II  It  II II II II  II  II  II  II  It  n  It  It  II  ti  n  n  n  H  n  n  n  ii  ii  ii  ii  ii  ii  ii  n  n  ti  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  n  n  n 

SUBROUTINE  fftfilter  (signal,  npts,  samhz,  corner,  order) 

11  If  II  II  II  II  II  II  II  II  n  II  II  II  it  II  II  n  11  n  ii  ii  it  ii  ii  ti  ii  ii  n  n  it  ii  ii  it  n  n  it  it  ti  n  n  n  ti  n  n  tt  ii  n  ii  ii  it  ti  ii  n  it  ii  n  ii  ii  ii  n  it  ii  ii  it  ii  ii  ii  ii  ii  ii  n 

"  This  subroutine  designs  a  Butterworth  filter  and  applies  it  to  the 
"  the  signal.  It  requires  a  fast  Fourier  transform  (FFTO  routine  (not 
"  provided  here)  because  filtering  is  done  in  the  frequency  domain. 

II 

"  Parameters ; 

II 

"  signal  . . .  array  containing  signal  before  and  after  filtering 

"  npts  . . .  number  of  samples  in  signal 

"  samhz  . . .  sampling  rate  of  signal,  in  Hertz 

"  corner  ...  -3  dB  corner  of  desired  Butterworth  filter,  in  Hertz 

"  npoles  . . .  number  of  poles  of  Butterorth  filter,  must  be  even 

II 

II  II  II  II  II  II II II  It  II  II  II II II  II  II  II  II  II  II  II  II  II  II  II  II II  It  II II II II II II  II  II  n  II 11  tt  II  n  n  n  It  ti  ii  ii  it  ii  n  n  ii  n  ii  ii  ii  it  ii  it  ii  ii  ii  ii  ti  n  n  n  ii  n  n  h 

REAL* 4  signal (*) 

funhz  =  samhz  /  npts 
wfund  =  funhz  /  corn 
order  =  nploes  /  2 
power  =  2  *  order 
nfreq  =  npts  /  2 

CALL  fft  (npts,  signal,  +1)  1  transform  to  frequency  domain 

DO  k  =  1,  nfreq 
ib  =  2  *  k 
ia  =  ib  -  1 

wfrq  =  (  k-1  )  *  wfund 

ginv  =  1  +  wfrq  **  power 

gain  =  1  /  SQRT (  ginv  ) 

signal (  ia  )  =  gain  *  signal (  ia  ) 

signal (  ib  )  =  gain  *  signal (  ib  ) 

END  DO 

CALL  fft  (npts,  signal,  -1)  I  transform  back  to  time  domain 

RETURN 

END 
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Program  to  design  Butterworth  low-pass  filters. 


n  If  n  II  II  II  II  II  II  n  ii  ii  ii  ii  ii  ii  n  n  ii  ii  ii  ii  ii  ii  ii  ii  ii  it  n  ii  ii  ii  ii  ii  ii  ii  ii  n  n  n  it  tt  n  ii  ii  ii  ii  ti  n  ii  n  ii  ii  ii  ii  ii  ii  n  ii  ii  ti  ii  ii  ii  ii  ii  ti  ii  ii  ii  ii  ii  ii  ii  ii  ii  n 

SUBROUTINE  designbutter  (samhz,  corner,  nsect,  acof,  bcof) 

It  II  ti  II II II  11  n  II  ti  II  II  It  It  ft  II II II 11 II II  II II II 11 II 11  ti  II  II  II  II  It  11 II II II II II  It  n  If  It  It  II II  II  II  II II II II II  It  II II II II II II 11  II  II  It  II II II  II  II II II 11  It  11 II II  II 
II 

"  Subroutine  to  design  low-pass  Butterworth  digital  filters.  The  filter  is 
"  obtained  by  using  the  bilinear  transformation  to  transform  analog  filter 
"  equations  to  digital  domain.  Filtering  is  accomplished  by  a  cascade  of 
"  second-order  sections  which  are  defined  by  the  order  of  the  filter. 

"  Implementation  in  the  time-domain  is  recursive.  Arguments  are: 

II 

"  Input : 

II 

"  samhz  ...  given  sampling  rate  (Hz)  of  digital  signal. 

II 

"  corner  . . .  given  filter  corner  frequency  (Hz)  where  the  magnitude 

"  is  -3  dB  (half-power  point) . 

II 

"  nsect  . . .  given  number  of  2nd-order  sections  (pole-pairs) .  The 

"  nxamber  of  poles  of  the  filter  will  be  2  x  nsect. 

tt 

"  Output : 

If 

"  acof  ...  coefficients  (A0,A1,A2)  of  2nd-order  filter  sections 

"  bcof  ...  coefficients  (B0,B1,B2)  of  2nd-order  filter  sections 

n 

"  Implementation: 

II 

"  Recursive  filtering  through  each  2nd-order  section  is  performed  by 

"  the  difference  ec[uation: 

11 

'■  y(n)  =  AO  *  X(n)  +  A1  *  X(n-l)  +  A2  *  X(n-2)  -  B1  *  Y(n-l)  -  B2  *  y(n-2) 

It 

II  II  II  11  II  II  II  II  II  II  n  II  11  II  II  It  II  II  n  II  II  II  n  n  ti  n  If  II  II  II  n  11  n  n  n  If  If  n  n  tt  II  It  n  n  n  n  II  n  n  n  n  It  It  II  II  II  II  n  II  n  n  It  II  II  ii  ii  ii  ii  ii  ii  it  n  ii  ii  ii  ii 

REAL*4  acof  (3,*),  bcof (3,*) 

REAL*4  pie  /3 . 1415926535/ 

wc  =  corner  /  samhz 
fact  =  TAN {  pie  *  wc  ) 
npoles  =  2  *  nsect 
sector  =  pie  /  npoles 
wedge  =  sector  /  2 . 


23 


DO  m  =  1,  nsect 


ang  =  wedge  *  (  2*m  -  1  ) 

xm  =  -  fact  *  COS {  ang  ) 
ym  =  fact  *  SIN(  ang  ) 

den  =  (  1.  -  xm  )**2  +  ym**2 

urn  =  (1.  -  xm**2  -  ym**2  )/  den 
vm  =  (  2 .  *  ym  ) /  den 

bcof(l,m)  =  1. 
bcof(2,m)  =  -2.  *  urn 
bcof(3,m)  =  um  *  urn  +  vm  *  vm 

sum  =  bcof(l,m)  +  bcof(2,m)  +  bcof{3,m) 

acof(l,m)  =  sum  /  4. 

acof(2,m)  =  stun  /  2. 

acof(3,m)  =  sum  /  4, 

END  DO 

RETURN 

END 
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Examples  of  filters  designed  with  the  design  butter  routine  listed  in  this  appendix 


Design  method:  Standard  Butterworth  Filter  corner  =  100  Hz 


Sampling 

rate :  samhz 

=  10000  Hz 

No.  sections  =  2 

Filter 

coefficients : 

Section  1 

Section  2 

aO  .  . 

acof(l)  = 

,932544E-03 

.963479E-03 

al  .  . 

• 

acof{2)  = 

.186509E-02 

.192696E-02 

3.2  •  • 

• 

acof(3)  = 

.932544E-03 

.963479E-03 

bl  .  . 

• 

bcof(2)  = 

-1.88661 

-1.94922 

b2  .  . 

bcof(3)  = 

.890340 

.953070 

Design 

method:  Standard 

Butterworth 

Filter  corner  =  300  Hz 

Sampling 

rate:  samhz  = 

10000  Hz 

No.  sections  =  2 

Filter 

coefficients : 

Section  1 

Section  2 

aO  .  . 

aco£(l)  = 

.754943E-02 

.826380E-02 

al  .  . 

aco£  (2 )  = 

.150989E-01 

.165276E-01 

a2  .  . 

aco£ (3)  = 

.754943E-02 

.826380E-02 

bl  .  . 

bco£(2)  = 

-1.67466 

-1.83313 

b2  .  . 

bco£(3)  = 

.704859 

.866181 

Design 

method:  Standard 

Butterworth 

Filter  corner  =  1000  Hz 

Sampling 

rate :  samhz  = 

10000  Hz 

No.  sections  =  2 

Filter 

coefficients : 

Section  1 

Section  2 

..  aO  .  . 

• 

aco£(l)  = 

.618852E-01 

.779563E-01 

al  .  . 

• 

aco£ (2 )  = 

.123770 

.155913 

a2  .  . 

• 

aco£(3)  = 

.  618852E-01 

.779563E-01 

bl  .  . 

• 

bco£(2)  = 

-1.04860 

-1.32091 

b2  .  . 

bcof(3) 

.296140 

.632739 

Design 

method:  Standard 

Butterworth 

Filter  corner  =  1650  Hz 

Sampling 

rate :  samhz  = 

16500  Hz 

No.  sections  =  2 

Filter 

coefficients : 

Section  1 

Section  2 

aO  .  . 

aco£(l)  = 

.618852E-01 

.779563E-01 

al  .  . 

aco£(2)  = 

.123770 

.155913 

a2  .  . 

aco£ (3 )  = 

.618852E-01 

.779563E-01 

bl  .  . 

bco£(2)  = 

-1.04860 

-1.32091 

b2  .  . 

bco£{3)  = 

.296140 

.632739 
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Appendix  C. 


Program  to  implement  a  second-order  filter  in  the  time  domain. 


II  II  II  II  II  II  II  II  II  II  ti  II  II  II  11  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  It  II  M  n  It  II  n  II  It  n  II  II  It  ii  ii  ii  it  ti  ti  it  ii  n  h  n  n  n  n  n  n  u  n  n  h  h  n  n  n  n  n  n  n  j|  n 

SUBROUTINE  f ilter_2nd_order  (x,  npt,  a,  b) 

11  II  II  It  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  11  11  II  11  II  II  II  II  II  II  II  II  II  II  II  II  If  If  If  f|  f|  II  If  f|  II  II  II  If  If  11  II  II  If  ,1  ,1  I,  „  II  I,  II  II  II  II  II  II  II  II  II  II  II  II  II  II 
II 

"  Subroutine  for  recursive  application  of  second-order  filter  to  a  time 
"  domain  signal . 

11 

"  Inside  this  routine,  filtering  is  forward.  Backward  filtering  may  be 
"  accomplished  by  reversing  the  signal  prior  to  calling  this  routine, 

"  then  restoring  the  order  upon  return  of  the  filtered  signal. 

II 

"  The  first  two  points  are  reflected  about  the  initial  point  to  produce 
''  a  reasonable  starting  method.  Other  initial  conditions  may  dictate 
"  other  starting  methods. 

II 

"  By  sliding  the  filter  window  along  the  time  axis,  the  need  for  auxiliary 
"  storage  is  eliminated,  allowing  the  full  utilization  of  computer  memory. 

II 

"  This  routine  illustrates  the  correct  usage  of  the  filter  coeffients 
"  in  the  difference  equation: 

II 

"  Y(n)  =  A0*X(n)  +  Al*X(n-l)  +  A2*X(n-2)  -  Bl*y(n-1)  -  B2*y(n-2) 

n 


"  Arguments : 

n  * 

"  xO 

II 

II 

"  npt 

It 

"  a{) 

"  b() 

II 

n 


upon  entry,  an  array  containing  the  unfiltered  signal, 
and  replaced  by  the  filtered  signal  upon  return. 

number  of  samples  in  the  x()  signal  array. 

array  containing  AO,  Al,  and  A2  coefficients  of  filter 
array  containing  BO,  Bl,  and  B2  coefficinets  of  filter 
Note:  BO  must  be  supplied  even  through  not  used. 


II  II  II  II  II  II  n  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  n  II  II  II  II  n  II  II  n  n  n  n  n  II  n  n  n  n  n  II  ii  ii  ii  ii  ii  ii  n  n  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii  ii 


REAL*4  x(*) ,  a(*) ,  b{*) 

aO  =  a(l) 
al  =  a(2} 
a2  =  a(3) 
bl  =  b(2) 
b2  =  b(3) 
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xnO  =  x(l) 

xnl  =  2  *  xnO  -  x(2) 
xn2  =  2  *  xnO  -  x(3) 

ynl  =  xn2 
yn2  =  xnl 

ynO  =  aO  *  xnO  +  al  *  xnl  +  a2  *  xn2  -  bl  *  ynl  -  b2  *  yn2 
x(l)  =  ynO 

xn2  =  xnl 
xnl  =  xnO 
xnO  =  x(2) 

yn2  =  ynl 
ynl  =  ynO 

ynO  =  aO  *  xnO  +  al  *  xnl  +  a2  *  xn2  -  bl  *  ynl  -  b2  *  yn2 
x(2)  =  ynO 

ynl  =  x(2) 
yn2  =  x(l) 

DO  n  =  3 ,  npt 

yn  =  aO  *  x(n)  +  al  *  x(n-l)  +  a2  *  x{n-2)  -  bl  *  ynl  -  b2  *  yn2 

x(n-2)  =  yn2 
yn2  =  ynl 
ynl  =  yn 

END  DO 

RETURN 

END 
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Appendix  D. 


Program  to  design  a  filter  per  SAE  J211  (draft)  guidelines. 


If  II  II  II  II  II  II  II  II  II  II  II  II  II  II  If  II  II  II  II  II  II  II  II  II  If  II  If  II  II  II  II  ti  II  II  It  It  n  It  It  ft  ti  ti  II  II  ti  II  ti  If  If  ti  It  11  II  II  II  II  II  If  II  II  It  II  II  II  II  II  II  II  II  If  11  II  It  II  II 

SUBROUTINE  design_J211  (samhz,  corner,  nsect,  aco£,  bcof) 

II  II  11  If  It  If  It  II  II  II  II  II  II  II  II  II  II  II  II  II  II  It  II  II  II  II  II  II  II  If  II  II  If  11  II  II  II  It  n  H  It  II II  II  It  II  II  II  It  11  If  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II 
II 

"  Subroutine  to  design  low-pass  digital  filter  using  equations  recommended 
"  by  the  SAE  J211  guideline (1994 ,  draft).  All  J211  channel  class  filters 
"  (CFG  60,  180,  300,  and  1000)  are  derived  from  a  4-th  order  Butterworth 
"  filter,  modified  to  keep  the  filter  response  inside  a  desired  frequency 
"  response  corridor.  Arguments  are: 

II 


Input : 
samhz 

corner 

nsect 


Output : 

acof 

bcof 


given  sampling  rate  (Hz)  of  digital  signal. 

. .  given  filter  corner  frequency  (Hz)  where  the  magnitude 
is  -3  dB  (half-power  point) .  This  is  equal  to  the  class 
CFG  of  the  filter,  divided  by  0.6  factor. 

.  .  given  ntimber  of  2nd-order  sections  (pole-pairs)  .  For  the 
J211  filters,  there  are  2  identical  sections. 


coefficients  (A0,A1,A2)  of  2nd-order  filter  sections 
...  coefficients  (B0,B1,B2)  of  2nd-order  filter  sections 


"  Implementation: 

n 

"  Recursive  filtering  through  each  2nd-order  section  is  performed  by 

"  the  difference  equation: 

"  Y(n)  =  AO  *  X(n)  +  A1  *  X(n-l)  +  A2  *  X(n-2)  -  B1  *  y(n-l)  -  B2  *  y(n-2) 


n  II  II  II  II  II  II  II  II  II  II  II  II  It  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  If  II  II  n  n  n  n  n  n  II II II II  n  n  n  n  ii  ii  ii  ii  ii  ii  n  n  n  n  n  n  n  n  n  n  h  h  n  n  n  n  h  h  n  n  n 


REAL*4  acof (3,*),  bcof (3,*) 

REAL*4  pie  /3 . 1415926535/ 
class  =  0.6  *  corner 

IF  (  NINT  (corner)  .EQ.  1650  )  class  =  1000 
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ts  =  1.00  /  samhz 

sr2  =  SQRT  (2.0) 

wd  =  2.d0  *  pie  *  class  *  2. 

arg  =  wd  *  ts  /  2 . 
wa  =  TAN  (arg) 
wa2  =  wa  *  wa 

den  =  (  1 .  +  sr2  *  wa  +  wa2 

aO  =  wa2  /  den 
al  =  2 •  *  aO 
a2  =  aO 

bO  =  1.0 

bl  =  2 .  *  (  wa2  -  1 .  )  /  den 
b2  =  (  1 .  -  sr2  *  wa  +  wa2  ) 


DO  m  =  1,  nsect 

aco£  (l/in) 

=  aO 

acof (2 ,m) 

=  al 

acof (3 ,m) 

=  a2 

bco£  (l/in) 

=  1. 

bco£ (2 ,m) 

=  bl 

bco£ (3 ,m) 

=  b2 

END  DO 

RETURN 

END 

0775 


) 


/  den 


Examples  of  filters  designed  with  the  designj211  routine  listed  in  this  appendix 


Design 

method:  SAE  J211 

guideline 

Sampling 

rate :  samhz  = 

10000  Hz 

Filter 

coefficients : 

Section  1 

aO  .  .  . 

aco£(l)  = 

.145237E-02 

al  .  . 

• 

acof(2)  = 

.290474E-02 

a2  .  . 

. 

acof (3 )  = 

.145237E-02 

bl  .  . 

bcof (2)  = 

-1.88934 

b2  .  . 

• 

bcof(3)  = 

.895153 

Design  method;  SAE  J211 
Sampling  rate:  samhz  = 

Filter  coefficients: 

guideline 
10000  Hz 

Section  1 

aO  .  . 

acof(l)  = 

.117963E-01 

al  .  . 

acof  (2)  = 

.235925E-01 

a2  •  • 

acof (3 )  = 

.117963E-01 

bl  .  . 

bcof (2)  = 

-1.67012 

b2  .  . 

bcof  (3)  = 

.717303 

Design  method:  SAE  J211 
Sampling  rate:  samhz  = 

Filter  coefficients: 

guideline 
10000  Hz 

Section  1 

aO  .  . 

• 

acof(l)  = 

.971846E-01 

al  .  . 

. 

acof  (2)  = 

.194369 

a2  .  . 

• 

acof (3)  = 

.971846E-01 

bl  .  . 

. 

bcof  (2 )  = 

- .945574 

b2  .  . 

• 

bcof(3)  = 

.334312 

Design  method:  SAE  J211 
Sampling  rate:  samhz  = 

Filter  coefficients: 

guideline 
16500  Hz 

Section  1 

aO  .  . 

a 

acof  (1)  = 

.987937E-01 

al  .  . 

a 

acof (2)  = 

.197587 

a2  .  a 

a 

acof (3)  = 

.987937E-01 

bl  a  a 

a 

bcof  (2 )  = 

- .935632 

b2  a  a 

a 

bcof  (3)  = 

.330807 

Filter  corner  =  100  Hz 
No.  sections  =  2 

Section  2 

.145237E-02 

.290474E-02 

.145237E-02 

-1.88934 

.895153 


Filter  corner  =  300  Hz 
No.  sections  =  2 

Section  2 

.117963E-01 

.235925E-01 

.117963E-01 

-1.67012 

.717303 


Filter  corner  =  1000  Hz 
No.  sections  =  2 

Section  2 

.971846E-01 

.194369 

.971846E-01 

-.945574 

.334312 


Filter  corner  =  1650  Hz 
No.  sections  =  2 

Section  2 

.987937E-01 

.197587 

.987937E-01 

-.935632 

.330807 
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Initial  distribution 


Commander,  U.S.  Army  Natick  Research, 
Development  and  Engineering  Center 
ATTN:  SATNC-MIL  (Documents 
Librarian) 

Natick,  MA  01760-5040 
Chairman 

National  Transportation  Safety  Board 
800  Independence  Avenue,  S.W. 
Washington,  DC  20594 

Commander 

10th  Medical  Laboratory 
ATTN:  Audiologist 
APO  New  York  09180 

Naval  Air  Development  Center 
Technical  Information  Division 
Technical  Support  Detachment 
Warminster,  PA  18974 

Commanding  Officer,  Naval  Medical 
Research  and  Development  Command 
National  Naval  Medical  Center 
Bethesda,  MD  20814-5044 

Deputy  Director,  Defense  Research 
and  Engineering 
ATTN:  Military  Assistant 
for  Medical  and  Life  Sciences 
Washington,  DC  20301-3080 

Commander,  U.S.  Army  Research 
Institute  of  Environmental  Medicine 
Natick,  MA  01760 

Library 

Naval  Submarine  Medical  Research  Lab 
Box  900,  Naval  Sub  Base 
Groton,  CT  06349-5900 


Executive  Director,  U.S.  Army  Human 
Research  and  Engineering  Directorate 
ATTN:  Technical  Library 
Aberdeen  Proving  Ground,  MD  21005 

Commander 

Man-Machine  Integration  System 
Code  602 

Naval  Air  Development  Center 
Warminster,  PA  18974 

Commander 

Naval  Air  Development  Center 
ATTN:  Code  602-B 
Warminster,  PA  18974 

Commanding  Officer 

Armstrong  Laboratory 

Wright-Patterson 

Air  Force  Base,  OH  45433-6573 

Director 

Army  Audiology  and  Speech  Center 
Walter  Reed  Army  Medical  Center 
Washington,  DC  20307-5001 

Commander/Director 
U.S.  Army  Combat  Surveillance 
and  Target  Acquisition  Lab 
ATTN:  SFAE-IEW-JS 
Fort  Monmouth,  NJ  07703-5305 

Director 

Federal  Aviation  Administration 
FAA  Technical  Center 
Atlantic  City,  NJ  08405 

Director 

Walter  Reed  Army  Institute  of  Research 
Washington,  DC  20307-5100 
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Commander,  U.S.  Army  Test 
and  Evaluation  Command 
Directorate  for  Test  and  Evaluation 
ATTN:  AMSTE-TA-M  (Human  Factors 
Group) 

Aberdeen  Proving  Ground, 

MD  21005-5055 

Naval  Air  Systems  Command 
Technical  ^r  Library  95  OD 
Room  278,  Jefferson  Plaza  II 
Department  of  the  Navy 
Washington,  DC  20361 

Director 

U.S.  Army  Ballistic 
Research  Laboratory 
ATTN;  DRXBR-OD-ST  Tech  Reports 
Aberdeen  Proving  Ground,  MD  21005 

Commander 

U.S.  Army  Medical  Research 
Institute  of  Chemical  Defense 
ATTN;  SGRD-UV-AO 
Aberdeen  Proving  Ground, 

MD  21010-5425 

Commander 

USAMRMC 

ATTN:  SGRD-RMS 

Fort  Detrick,  Frederick,  MD  21702-5012 

HQ  DA  (DASG-PSP-O) 

5109  Leesburg  Pike 
Falls  Church,  VA  22041-3258 

Harry  Diamond  Laboratories 
ATTN:  Technical  Information  Branch 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 


U.S.  Army  Materiel  Systems 
Analysis  Agency 

ATTN:  AMXSY-PA  (Reports  Processing) 
Aberdeen  Proving  Ground 
MD  21005-5071 

U.S.  Army  Ordnance  Center 
and  School  Library 
Simpson  Hall,  Building  3071 
Aberdeen  Proving  Ground,  MD  21005 

U.S.  Army  Environmental 
Hygiene  Agency 
ATTN:  HSHB-MO-A 
Aberdeen  Proving  Ground,  MD  21010 

Technical  Library  Chemical  Research 
and  Development  Center 
Aberdeen  Proving  Ground,  MD 
21010-5423 

Commander 

U.S.  Army  Medical  Research 
Institute  of  Infectious  Disease 
ATTN:  SGRD-UIZ-C 
Fort  Detrick,  Frederick,  MD  21702 

Director,  Biological 
Sciences  Division 
Office  of  Naval  Research 
600  North  Quincy  Street 
Arlington,  VA  22217 

Commandant 
U.S.  Army  Aviation 
Logistics  School  ATTN:  ATSQ-TDN 
Fort  Eustis,  VA  23604 

Headquarters  (ATMD) 

U.S.  Army  Training 
and  Doctrine  Command 
ATTN:  ATBO-M 
Fort  Monroe,  VA  23651 
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lAF  Liaison  Officer  for  Safety 
USAF  Safety  Agency/SEFF 
9750  Avenue  G,  SE 
Kirtland  Air  Force  Base 
NM  87117-5671 

Naval  Aerospace  Medical 
Institute  Library 
Building  1953,  Code  03L 
Pensacola,  FL  32508-5600 

Command  Surgeon 
HQ  USCENTCOM  (CCSG) 

U.S.  Central  Command 
MacDill  Air  Force  Base,  FL  33608 

Director 

Directorate  of  Combat  Developments 

ATTN;  ATZQ-CD 

Building  515 

Fort  Rucker,  AL  36362 

U.S.  Air  Force  Institute 
of  Technology  (AFIT/LDEE) 
Building  640,  Area  B 
Wright-Patterson 
Air  Force  Base,  OH  45433 

Henry  L.  Taylor 
Director,  Institute  of  Aviation 
University  of  Illinois-Willard  Airport 
Savoy,  IL  61874 

Chief,  National  Guard  Bureau 
ATTN:  NGB-ARS 
Arlington  Hall  Station 
111  South  George  Mason  Drive 
Arlington,  VA  22204-1382 

AAMRL/HEX 

Wright-Patterson 

Air  Force  Base,  OH  45433 


Commander 

U.S.  Army  Aviation  and  Troop  Command 
ATTN:  AMSAT-R-ES 
4300  Goodfellow  Bouvelard 
St.  Louis,  MO  63120-1798 

U.S.  Army  Aviation  and  Troop  Command 
Library  and  Information  Center  Branch 
ATTN:  AMSAV-DIL 
4300  Goodfellow  Boulevard 
St.  Louis,  MO  63120 

Federal  Aviation  Administration 
Civil  Aeromedical  Institute 
Library  AAM-400A 
P.O.  Box  25082 
Oklahoma  City,  OK  73125 

Commander 

U.S.  Army  Medical  Department 
and  School 
ATTN:  Library 
Fort  Sam  Houston,  TX  78234 

Commander 

U.S.  Army  Institute  of  Surgical  Research 

ATTN:  SGRD-USM 

Fort  Sam  Houston,  TX  78234-6200 


Air  University  Library 
(AUL/LSE) 

Maxwell  Air  Force  Base,  AL  36112 

Product  Manager 

Aviation  Life  Support  Equipment 

ATTN:  SFAE-AV-LSE 

4300  Goodfellow  Boulevard 

St.  Louis,  MO  63120-1798 
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Commander  and  Director 

USAE  Waterways  Experiment  Station 

ATTN:  CEWES-IM-MI-R, 

CD  Department 
3909  Halls  Ferry  Road 
Vicksburg,  MS  39180-6199 

Commanding  Officer 
Naval  Biodynamics  Laboratory 
P.O.  Box  24907 
New  Orleans,  LA  70189-0407 

Assistant  Commandant 
U.S.  Army  Field  Artillery  School 
ATTN:  Morris  Swott  Technical  Library 
Fort  Sill,  OK  73503-0312 

Mr.  Peter  Seib 

Human  Engineering  Crew  Station 
Box  266 

Westland  Helicopters  Limited 
Yeovil,  Somerset  BA20  2YB  UK 

U.S.  Army  Dugway  Proving  Ground 
Technical  Library,  Building  5330 
Dugway,  UT  84022 

U.S.  Army  Yuma  Proving  Ground 
Technical  Library 
Yuma,  AZ  85364 

AFFTC  Technical  Library 
6510  TW/TSTL 
Edwards  Air  Force  Base, 

CA  93523-5000 

Commander 
Code  3431 

Naval  Weapons  Center 
China  Lake,  CA  93555 


Aeromechanics  Laboratory 
U.S.  Army  Research  and  Technical  Labs 
Ames  Research  Center,  M/S  215-1 
Moffett  Field,  CA  94035 

Sixth  U.S.  Army 
ATTN:  SMA 

Presidio  of  San  Francisco,  CA  94129 
Commander 

U.S.  Army  Aeromedical  Center 
Fort  Rucker,  AL  36362 

Strughold  Aeromedical  Library 

Document  Service  Section 

2511  Kennedy  Circle 

Brooks  Air  Force  Base,  TX  78235-5122 

Dr.  Diane  Damos 
Department  of  Human  Factors 
ISSM,  use 

Los  Angeles,  CA  90089-0021 

U.S.  Army  White  Sands 
Missile  Range 
ATTN:  STEWS-IM-ST 
White  Sands  Missile  Range,  NM  88002 

Director,  Airworthiness  Qualification  Test 
Directorate  (ATTC) 

ATTN:  STEAT-AQ-O-TR  (Tech  Lib) 

75  North  Flightline  Road 

Edwards  Air  Force  Base,  CA  93523-6100 

Ms.  Sandra  G.  Hart 
Ames  Research  Center 
MS  262-3 

Moffett  Field,  CA  94035 

Commander 

USAMRMC 

ATTN:  SGRD-UMZ 

Fort  Detrick,  Frederick,  MD  21702-5009 
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Commander 

U.S.  Army  Health  Services  Command 

ATTN:  HSOP-SO 

Fort  Sam  Houston,  TX  78234-6000 

U.  S.  Army  Research  Institute 
Aviation  R&D  Activity 
ATTN;  PERI-IR 
Fort  Rucker,  AL  36362 

Commander 

U.S.  Army  Safety  Center 
Fort  Rucker,  AL  36362 

U.S.  Army  Aircraft  Development 
Test  Activity 
ATTN:  STEBG-MP-P 
Cairns  Army  Air  Field 
Fort  Rucker,  AL  36362 

Commander 

USAMRMC 

ATTN:  SGRD-PLC  (COL  R.  Gifford) 
Fort  Detrick,  Frederick,  MD  21702 

TRADOC  Aviation  LO 
Unit  21551,  Box  A-209-A 
APO  AE  09777 

Netherlands  Army  Liaison  Office 

Building  602 

Fort  Rucker,  AL  36362 

British  Army  Liaison  Office 

Building  602 

Fort  Rucker,  AL  36362 

Italian  Army  Liaison  Office 

Building  602 

Fort  Rucker,  AL  36362 


Directorate  of  Training  Development 

Building  502 

Fort  Rucker,  AL  36362 

Chief 

USAHEL/USAAVNC  Field  Office 
P.  O.  Box  716 

Fort  Rucker,  AL  36362-5349 

Commander,  U.S.  Army  Aviation  Center 
and  Fort  Rucker 
ATTN:  ATZQ-CG 
Fort  Rucker,  AL  36362 

Dr.  Sehchang  Hah 
Dept,  of  Behavior  Sciences  and 
Leadership,  Building  601,  Room  281 
U.  S.  Military  Academy 
West  Point,  NY  10996-1784 

Canadian  Army  Liaison  Office 

Building  602 

Fort  Rucker,  AL  36362 

German  Army  Liaison  Office 

Building  602 

Fort  Rucker,  AL  36362 

French  Army  Liaison  Office 
USAAVNC  (Building  602) 

Fort  Rucker,  AL  36362-5021 

Australian  Army  Liaison  Office 
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